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Abstract 

The entanglement entropy in a quantum field theory between two regions of space has been shown 
in simple cases to be proportional to the volume of the hypersurface separating the regions. We 
prove that this is true for a free scalar field in an arbitrary geometry with purely spatial curvature 
and obtain a complete asymptotic expansion for the entropy. 
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I. INTRODUCTION 



Modern developments in efforts to consistently combine gravity and quantum mechanics 
have indicated that quantum field theory has too many degrees of freedom. The entropy 
of a geometrical object seems to depend on the area of a boundary surface. In the case 
of a black hole, thermodynamic arguments suggest that its entropy is proportional to the 
area of its event horizon. This expectation has been confirmed on the microscopic side by 
calculations for some special cases of black holes. 

The separation of a system into two subsystems gives rise to the notion of the entangle- 
ment entropy which quantifies the quantum correlations between the two subsystems. One 
can take the subsystems to be two regions of space separated by a boundary surface. Early 
calculations in quantum field theories indicated that the entanglement entropy between the 
degrees of freedom in two regions separated by a boundary surface is proportional to the 
area of the boundary surface , |4|. These developments have led to the intuition 

that the entanglement entropy is dominated by degrees of freedom close to the boundary 
surface, so it is natural to expect that the entanglement entropy is proportional to the area. 
The connection between the entropy of a geometric system, on the one hand, and the en- 
tanglement entropy between the quantum field theory degrees of freedom between spatially 
disconnected regions, on the other hand, is not clear. There have even been suggestions that 
the black hole entropy is entirely entanglement entropy [5]. This suggests the possibility 
that when the boundary surface is taken to be an event horizon, the two types of entropy 
are identical or at least related. Therefore it becomes important to examine the properties 
of entanglement entropy in more general cases and attempt to understand its properties 
when the entanglement involves regions separated by a horizon. 

The entanglement entropy is ultraviolet divergent and must be regularized. Presumably 
a more fundamental theory at the Planck scale will provide the mechanism that eliminates 
the divergence. Since we are not yet aware of how the underlying theory regularizes the 
divergence, we are forced to do so by hand in calculations. The regularization involves 
the Planck scale and the expectation that entanglement entropy may be connected to the 
theory of quantum gravity is suggested by simple dimensional analysis. If indeed the entan- 
glement entropy is proportional to the area of the boundary surface in Planck units, then 
the connection to a holographic principle similar to the one suggested by the AdS/CFT 
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correspondence {(], Q, 0, Q , may apply to the entanglement entropy 

After the initial studies of the entanglement entropy , most research in the 

subject has been devoted to CFTs. When a field theory is conformal, there are additional 
tools to compute the entanglement entropy, which, in fact, is proportional to the central 
charge 11]. Even when the CFT is modified by adding mass deformations 11, Q, 
the same property holds. Also, recent calculations of entanglement entropy in this context 
suggest a holographic interpretation Q, [ijj]. The connection between the entanglement 



entropy and holography has also been discussed in more general contexts [16l . |17] . Scaling of 
the entablement entropy with the area of the boundary have also been verified by numerical 

nn 

computations [18|, [19[ . 

In this paper, we investigate the entanglement entropy of a quantum field theory in the 
case of an arbitrary boundary surface embedded in a background with spatial curvature. 
Our results contain many previously derived results as special cases, but are more general 
because we consider arbitrary geometries. Ultimately one would want to consider the even 
more general case involving spacetime curvature, so the comparison with other forms of 
entropy can be made explicit. 

II. ENTANGLEMENT ENTROPY 

In this section, we compute the entanglement entropy by using two main tools, the replica 
method and the heat kernel method. 

In the replica method (see, for example, Ref. 3|), the entropy is expressed as the limit 
k — > 1 of an expression involving the kth power of the (reduced) density matrix. Expressing 
the density matrix as a path integral, we are lead to consider the manifold which is the 
result of gluing k copies of the original manifold. This gives an expression of the entropy 
in terms of the path integral over closed curves in the glued manifold, which is described in 
the subsection III Al 

For a free scalar field on an arbitrary base manifold, the resulting path integral is expressed 
in terms of the spectral quantities of a differential operator. It is convenient to study these 
quantities by the heat kernel method. We use it in the subsection III Bl to obtain an expression 
for the entropy in which the dependence on the hypersurface is factored out. 

This leads to several properties of the entropy which we derive in the subsection III CI In 
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particular, we compute the entropy for a hypersurface which is a direct product of manifolds, 
prove the addititivity property of the entropy, and compute the leading terms of asymptotic 
expansions of the entropy. 



In the heat kernel method (see, for example, Refs. 20|, l2lj). the spectral information 
is obtained from an asymptotic expansion of the trace of the heat kernel of the operator. 
All terms in the expansion are determined by the geometry of the underlying manifold. 
Although they can be computed in principle, the computations are quite complicated in 
practice. The parameters in the resulting asymptotic expansion are an ultraviolet cutoff 
scale and geometric scales associated with the manifold. In the subsection III Dj. this leads 
to an asymptotic expansion for the entropy, which involves geometric quantities associated 
with the hypersurface. We show that the term proportional to the volume of the manifold 
is absent and the leading term is proportional to the volume of the hypersurface. 



A. Replica Method 

We consider a field theory on a generally curved space which is divided by an arbitrary 
hypersurface into two parts. The quantum fields in the two parts are entangled, and our goal 
is the calculation of the entanglement entropy. Let M be an n — 1 dimensional Riemannian 
manifold without a boundary and let £ C M be a closed hypersurface (a submanifold of 
codimension 1). £ divides M into two parts, the interior part M' and the exterior part 
M". Let be a field on M, and ((f)' , 0") its restrictions to (M',M"), and let be a 

wave function corresponding to the field having the value (0', 0") on (M', M"). The density 
matrix for (0', 0") is 

P (0' 1; 0'/, <p> 2 , 0' 2 ') = w v 0' 1 'M0 / 2 ,0' 2 / ), (i) 

and the reduced density matrix for 0' is obtained by tracing over the degrees of freedom of 
the field on M" , 




The reduced density matrix p' represents the mixed state with the associated entropy 

S' = - tr ( -^-log-^- ] - limf 1 - 4" ] logtrp /fc . (3) 
\trp' 6 trp7 k^i\ dk J 6 H w 
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The quantity tr p' = J d(f)' p' ((f)' , <$') in the denominator guarantees the correct normalization 
for the density matrix. The second e q oality emhodres the replica method (see, for exampie, 
Ret. i). 

In order to obtain the path integral representation for the fcth power of the density 
matrix we introduce an auxiliary field ip(r, x) defined on N = R x M and which satisfies the 
boundary condition ip(0,x) = <po(x). The parameter r represents Euclidean time. Let I(ip) 
be an action for the field <p. The wave function is ip(<f)) = Z(N, <fto, <ft), where 

Z(N,<f> ,<j))= [ dp exp (-J(p)) (4) 

JC{N,4>o,(t>) 

is a path integral over the space C(N,<j) ,<j)) of curves defined on iV and which satisfy 
boundary conditions cp(0,x) = <fio{x) and (p(T,x) = <fi(x) for some Tel. Using T = — oo 
for ip(4>[, 4>") and T = oo for ip(4>' 2 , <p"), we find 

pf(tiM = Z(N,<l>' 1 ,<l>' 2 ). (5) 

The function (p(r,x) has a discontinuity at r = since <p>(0~,x) = <j)[(x) and ip(0 + ,x) = 
4>' 2 {x). However, ip(r,x) is continuous on the manifold N±, which is defined as the manifold 
N with the cut along {r = 0} x M'. 
The kth power of the density matrix is 

p' fc (0;,0' fc+1 ) = J ^^•••# / fc p'(0 / 1 ,^)p / (0 / 2 ,0^)---p'^,0U 1 ). (6) 

Let R x M(i), ...,lx M(fc) be k copies of R x M. By cutting every R x along {r, = 
0} x M^s and gluing them in such a way that {fj = - } x AfL is identified with {r i+1 = 
+ } x M/' i+1 x for i = 1, . . . , k — 1, we obtain the manifold iV fc . See Fig. [H This gives 

pVi,0Ui) = ^X0Ui)- (7) 
Identifying {t^ = - } x M/ fc j with {t\ = + } x ML-,, we obtain the manifold A^^. This gives 

tr p' k = Z(N k )= [ dtpexp{-I(<p)), (8) 

JC(N k ) 

which is a path integral over all closed curves in iV^. This quantity gives the entanglement 
entropy for p' via Eq. ([3]) . If instead we were to trace the density matrix over the degrees of 
freedom in M', we would obtained the reduced density matrix p" for <fi". It is easy to show 
that the entanglement entropy for p" is the same, S' = S", and we denote the common value 
by Sx to emphasize its dependence on the surface E. 
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FIG. 1: The replica method involves cutting the original manifold N along {r = 0} x M' and 
gluing k such cut copies of N along {r = 0} x M' to form the manifold iV& with sheets. We 
identify = 0"} x ML with {r m = 0+} x M' (i+1) for i = 1, . . . , k - 1, and {r k = 0~} X M ( ' fc) 
with {n = + } x M/jq. When M = M and S is a point P, the construction gives the 2-dimensional 
cone manifold Ck = K + X S^, where S\ is the unit circle S 1 which is parametrized by < 9 < 2irk. 
The quantity 2ir{\ — k) is called the deficit angle. Note that Ck is the Riemann surface of the 
holomorphic function z i— > z k . 

B. Heat Kernel 

To proceed with an explicit computation, we choose the free scalar field with the action 

I(<p) = 2- 1 I co N cpD N cp, (9) 

J N 

where Dn = A^ + m 2 , A^ is the scalar Laplace operator for N, ujn is the volume form for 
N, and m is the mass of the field (p. Performing a Gaussian integral, we find 

Z(N k ) = Z^detD Nk y 1 / 2 , (10) 

where Z is a constant independent of k. Since D^ k is a non- negative elliptic operator, we 
can define its determinant by 

poo 

logdet D Nk — log det E Nk = — I dt t^ 1 (tr exp (—tD Nk ) — trexp (—tE Nk )), (11) 

Jo 

where E^ k is any other non- negative elliptic operator on N k - (To prove this equation, 
one writes the analogous equation relating eigenvalues of D^ k and E^ k .) The quantity 
exp (—tD Nk ) is called the heat kernel of the operator D Nk , and 

K(t,D Nk ) = trexp (-tD Nk ) (12) 

is its L 2 trace. We find K(t, D Nk ) = exp (—tm 2 )K(t, A Nk ). 



6 



The integral over t in Eq. (TTTT) diverges for small t. To obtain a finite result, we replace 
the lower limit of integration over t by a regularization parameter A 2 (an ultraviolet cutoff), 

POO 

/ dtr x K{t,D Nk ) =trr(0,A 2 AvJ. (13) 

J A 2 

Here T is the incomplete Gamma function which is given either by the integral representation 

roc 

F(a,z)= duu a ~ l ex.p(-u) (14) 

J z 

or by the series representation 



DC 



T(a,z) = T(a)-z a ^2 . 1 ^ a ^ 0, -1, -2, . . . , (15) 
z) = + 1) - lagz) -z~ l Y< ( ^ ., , J = 0,1,2,..., (16) 

where ^(/ + 1) = —7 + X]'=ii 1 an d 7 * s the Euler constant. We will later need the 
incomplete Gamma function for nonzero values of a as well. 

We choose E^ k to be a unit operator times a constant with the dimension of inverse 
length squared; this leads to vanishing of its contribution to the entropy. Similarly, the 
contribution from the constant Zq vanishes. The regularized entropy becomes 

S S (X) = 2- 1 lim(l - A) trT(0, X 2 D Nk ). (17) 

We can factor the dependence of Se(A) on £ proceeding as follows. Locally, N k = C*. x S, 
where Ck = K + x S\ is the 2-dimensional cone manifold, and S\ is the unit circle S 1 which 
is parametrized by < 9 < 2ixk. The quantity 27r(l — k) is called the deficit angle. Note 
that Ck is the Riemann surface of the holomorphic function z 1— > z k . Giving a product 
metric, we find 

A Nk = A Ck ® l s + l Ck ® A s , (18) 

which gives 

K(t,A Nk ) = K(t,A Ck )K(t,Ax). (19) 
This factorization reveals the special role played by the entropy for a point P, when M — R, 



5 P (A) = 2 



-1 



/"OO 

/ c/tr 1 exp(-tm 2 )C(t), (20) 
A 2 
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where 

c(«) = £s(i-!)*(*,a«). (21) 

A simple computation gives the expression 

Sx{\) = -f~d» A E ), (22) 

in which the dependence on E is factored out. This equation leads to several properties of 
the entropy, which we now derive. 



C. Several Properties of Entropy 

1. Let (r, 9) be local polar coordinates for C k , and £ > 0. Under the scaling transforma- 
tion (t, r, 9) i ► (£ 2 t, £r, 9), we have C k i-> C k , tA Ck h-> tA Ck , and so K(t, A C J h-> if (t, A C J. 
This implies C(t) = C — const and thus 

5p(A) =2~ 1 CT(0,A 2 m 2 ). (23) 

We will compute C in the next subsection and the appendix. 

2. Let (r = x 1 , x 2 ) be local coordinates for C k , let (x 3 , . . . , x n ) be local coordinates for 
E, and let £ > 0. Under the transformation A ^ £A, x J ' h- > j = 3, ...,n, we have 
A s i — > A 2i ^ = £~ 2 A S and <Se(A) i— ► 5s,^(£A), where 

s,AW--J !x & — (^n^J -flfpri" (24) 

Using Eq. ( |23l) . we find 

S S ,,(£A) = - jHdi/ exp (-(£ 2 - l)zAn 2 )^|^. (25) 
It follows that \im^ (dS E> ^X)/d^) = 0. Since vol (E) ^ £"~ 2 vol (E), this implies 

S E (A) ~C"A 2 - n vol(E), A^O, (26) 

where C" = const. 

3. Since is the smallest eigenvalue of As, we have K(t, As) ~ 1, t — > oo. This gives 

Ss(A) ~ 2- 1 CT(0, A 2 m 2 ), A -> oo. (27) 



S 



9(a- 1 Eiua- 1 s 2 ) 




FIG. 2: The hypersurfaces used in the formulation of the additivity property. 

Interestingly, this coincides with the expression in Eq. (|23j) for Sp(X) for arbitrary A. This 
can be understood as the result of the physical scales of £ becoming irrelevant as the cutoff 
A tends to infinity. 

4. Let Ei and E 2 be closed hypersurfaces in M. Let d and d~ l be operators defined by 
dM[ = Ej and <9 _1 £j = where M[ is a part of M inside E, for i — 1, 2. Being an integral 
over E, the quantity K(t, A E ) is linear in S. It follows that 

A Sl ) + K(t, A S2 ) = K(t, A a(a -i SlU9 -i S2) ) + A e(a -i Ein9 -i S2) ), (28) 

and thus the entropy satisfies the additivity property 

Sei(A) + Ss 2 (A) = 5 , 8 ( 9 -i SlUa -i E2 )(A) + <S' 9 ( 9 -i Ein9 -i Ea )(A). (29) 

See Fig. [2J For an arbitrary system, the entanglement entropy satisfies the strong subaddi- 
tivity property, which requires '>' instead of '=' in Eq. fl29|) . 



D. Asymptotics 

We now derive the asymptotic expansion of S^A) for A — > 0. It is clear that this requires 
knowledge of the asymptotic behavior of K(t, A E ) for t — > 0. For an n-dimensional manifold 
L, such an asymptotic is given by 

oo 

K(t,A L ) ~ £;^ l -"V2 0l ( Aj .) > t ^ o, (30) 

1=0 

where 

a t (A L ) = J uj L ai(x L ,A L ), (31) 
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and cli{xl, Al) are the heat kernel coefficients for Al- The above factorization of K(t, Ajy k ) 



leads to 



i 

ai(x Nk ,A Nk ) = ^2a j (x Ck ,A Ck )a l _ j (xx,A T: ). (32) 

3=0 



The coefficients ai(xL, Al) are completely determined by the geometry of L. For a mani- 
fold without boundary, ai(xi, Al) = for odd I. All coefficients ai(xL, Al) are polynomials 
in the covariant derivatives of the Riemann tensor (RL)abcd, the Ricci tensor (i?z) a 6, and the 
scalar curvature Rl of L. Explicit expressions for several first coefficients are available in 
the literature (see, for example, Ref. 20j). For example, 

a (* L ,A L ) = (47r)-"/ 2 , (33) 
a 2 (x L ,A L ) = (47r)- n / 2 6" 1 i? L , (34) 
a 4 (x L , Al) = (47r)-"/ 2 360- 1 {-12A L R L + 5R 2 L - 2 J2( R l) ab(R L ) ab 

a,b 

+ ^ ^ {RLlabcdiRhjabcd^j ■ (35) 



a,b,c,d 



It has been shown 



221 ] that the only nonzero heat kernel coefficients for Ck are 



a (x Cfc ,A c J = (47r)- 1 , (36) 
a 2 (x Ck ,A Ck ) = (47r)- 1 6- 1 47r(l - k)5 Ck , (37) 

where Sc h is the delta function at the origin of Ck- This gives C = 6 _1 . (In the appendix, 
we derive this result.) The regularized entropy becomes 

S P {\) = 12-^(0, A 2 m 2 ), (38) 
S S (X) = 12-Hr r(0,A 2 Ds3). (39) 

In terms of the integrated heat kernel coefficients of E, we find 

oo 

5 E (A) ~ 12- 1 ^m n -'- 2 r((2 + / - n)/2, (Am) 2 )a,(A s ), A -> 0. (40) 

1=0 

This asymptotic expansion is our main result. For A — * 0, Se(A) depends only on the spectral 
properties of the operator \ 2 D^. Equivalently, the entropy depends only on parameters m, 
A, and on geometric invariants associated with E. The asymptotic expansion for S's(A) 
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involves log Am and the powers of Am. The leading term in the entropy is 

Sp(X) ~ 12" 1 (-2 log Am -7), A -»• 0, (41) 
S s {\) ~ 12 _1 (n/2- l)- 1 A 2 -"(47r) 1 - n / 2 vol(S), A -»• 0. (42) 

The term of order A _n vol (N) in S's(A) is absent; it would be the extensive contribution to 
the entropy. 

We remark on the case n = 2. (See also Ref. [ljj] for a similar discussion.) The entangle- 
ment entropy for a critical 2-dimensional CFT with the central charge c is asymptotically 
S ~ (c/3) log (£/A), where I is the size of the system [ll]. For a massive theory with the 
correlation length £, the entropy becomes 5 ~ (c/6)z/ log (£/A) for £ ^> where z/ is the 
number of components of (zero dimensional) S. Setting c=l,i/ = l,£~ m -1 , we recover 
Eq. flUD. 



III. CONCLUSION 



We have calculated the asymptotic expansion of the entanglement entropy for a free scalar 
field in arbitrary background geometry. The expansion parameter is the ultraviolet cutoff A 
which is needed to regularize the entropy. We have found that the entropy depends only on 
geometric invariants associated with the boundary surface S. The extensive contribution to 
the entropy, the term of order A _n vol(iV), is absent. The leading term is proportional to 
A 2 - n vol(S). 

We have considered a situation with spatial curvature only and with time included in 
only a trivial way as a product. Our calculation does not utilize a spacetime that is a 
solution of Einstein's equations. Further research may involve extending this calculation to 
cases involving spacetime curvature. Another interesting direction to pursue is to include 
interactions since this, at least intuitively, can potentially change the area dependence of 
the entropy. The studies of interactions have been mostly limited to CFTs, but their role 
in QFT remains largely unexplored. These cases involve gravity more explicitly and may 
make a connection between holographic entropy and entanglement entropy more obvious. 
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APPENDIX 

Here we compute the heat kernel coefficients for the cone Let £ > 0. Under the trans- 
formation (t,r,6) i-> (£ 2 t,£r,#), we have K{t,A Ck ) ^ K{t,A Ck ), a/(A Cfc ) h-> £ 2 -'a*(A c J. 
Since Ck does not have a length scale associated with it, this implies ao(Ac k ) = oo, 
a2(Ac fc ) = const, a;(AcJ = 0, I > 4. a (xc k , Ac k ) is given by Eq. fl33l) . and to compute 
a 2 (xc t ) Aq.) from Eq. (1341) . we need to know the scalar curvature of Ck, with computation 
of which we now proceed. 

Ck is singular at r = if k ^ 1. We consider it as a limit = lim e ^ Cfc,e, where Cfe j£ 
is a regular manifold. On Cfc i£ we take an orthonormal frame (u; 1 ,^ 2 ) = (fdr,rd8), where 
the regularization function f(k,r,e) is an arbitrary smooth function satisfying conditions 

1 /2 

lim r ^o / — k, lim^o / = 1- An example of such a function is / = (fc 2 + (1 — /c 2 )(gr) £ ) , 
where £ > and q > is an arbitrary constant with the dimension of inverse length. In 
what follows, we proceed with arbitrary / satisfying the above conditions. 

Let lo and Q be the 2x2 antisymmetric matrices of connection and curvature forms. 
Cartan's equations for Cfc j£ , 



io 2 A rd6 



(A.l) 



dr Ad9 + uj 2 x A fdr 



0. 



(A.2) 



du> 2 



n 1 



2- 



(A.3) 



have the solution 



IO 2 



nde 



(A.4) 



rfgdr A c/6 1 , 



(A.5) 



where g = r 1 f ^{df /dr). The nonzero components of the Ricci tensor are (Rc ke )n 
(Rc k J22 = 9, and the scalar curvature is Rc k£ = 2g. 
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To obtain non-regularized quantities, we consider the limit e — > 0. Using lim r ^o / = k, 
lim e ^o / — 1) Wc t e — rfdr A for an arbitrary function h(r) satisfying h(oo) = 0, we find 

limy u CkE gh = 27rk\im(^-f- 1 h\ r r 2™+ dr f-\dh/dr)^j = 2tt(1 - fc)/i(0). (A.6) 

This implies limbos' — 27r(l — k)5c k , where 5c k is the delta function at the origin of Ck- 
Thus, Rc k = 47r(l — k)5c k , the only nonzero heat kernel coefficients for Ck are 

a (x Ck ,A Ck ) = (4n)- 1 , (A.7) 
a 2 (x Ck ,A Ck ) = (Any^-Hnil - k)5 Ck , (A.8) 

so that we may identify C = 6 _1 . Alternatively, since C(i) in Eq. (I2TI) is a constant, we can 
compute it by taking the limit t — > in the expansion 

1 -afe)E* (, " a)/SIa *( A ^)' ( A - 9 ) 

Since the coefficients a (xc fc , Ac fe ), ai(xc k , ApJ, Z > 4 do not contribute to C, we find 
C = 6" 1 . 
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